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ABSTRACT 

We introduce PIRATE, a new remotely-operable telescope facility for use in research and 
education, constructed from 'off-the-shelf hardware, operated by The Open University. We 
focus on the PIRATE Mark 1 operational phase where PIRATE was equipped with a widely- 
used 0.35m Schmidt-Cassegrain system (now replaced with a 0.425m corrected Dall Kirkham 
astrograph). Situated at the Observatori Astronomic de Mallorca, PIRATE is currently used to 
follow up potential transiting extrasolar planet candidates produced by the SuperWASP North 
experiment, as well as to hunt for novae in M31 and other nearby galaxies. It is operated by 
a mixture of commercially available software and proprietary software developed at the Open 
University. We discuss problems associated with performing precision time series photometry 
when using a German Equatorial Mount, investigating the overall performance of such 'off-the- 
shelf solutions in both research and teaching applications. We conclude that PIRATE is a 
cost-effective research facility, and also provides exciting prospects for undergraduate astronomy. 
PIRATE has broken new ground in offering practical astronomy education to distance-learning 
students in their own homes. 

Subject headings: Astronomical Instrumentation, Data Analysis and Techniques, Extrasolar Planets 



1. Introduction 

PIRATE (Physics Innovations Robotic Astro- 
nomical Telescope Explorer) is a remote telescope 
facility situated on the Balearic island of Mallorca 
(2°57'03.34"E, 39°38'34.31" N), at the Observa- 
tori Astronomic de Mallorca, ~162m above sea 
level. It is a remotely-operable facility used for 
research and undergraduate teaching, allowing for 
a practical astrophysics component to be imple- 
mented within The Open University's (OU, Mil- 
ton Keynes, UK) distance-learning modules. PI- 



RATE is constructed entirely from 'off-the-shelf 
hardware, and is operated primarily with commer- 
cial observatory control software, with some addi- 
tional proprietary software developed to support 
simultaneous use by groups of students. The total 
cost of the facility, combining all purchased hard- 
ware and software is of order $150k. 

The funding for PIRATE was primarily made 
available by a teaching innovation initiative 
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(piCETL \j) to explore the integration of a 
remotely-controlled observatory into university- 
level distance teaching modules. The aim of the 
project was to reproduce the hands-on experience 
of a traditional lab course with real-time access 
to a telescope from any computer linked to the 
internet. 

In this paper we refer to the PIRATE Mark 
1 facility (unless mentioned otherwise) which fea- 
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tured a Celestroii C14 0.35m Schmidt-Cassegrain 
Optical Tube Assembly (OTA). This OTA is cur- 
rently no longer in use within the PIRATE facil- 
ity, having been superseded in August 2010 by a 
0.425m PlaneWave System^ CDK17 f/6.8 Cor- 
rected Dall-Kirkham Astrograph telescope, which 
makes use of a new custom micro-focuser (pro- 
vided by PlaneWave), which we designate PI- 
RATE Mark 2. The Mark 1 hardware remains 
an obvious choice for cost-effective astronomy, so 
we report in detail on its capabilities, features, and 
drawbacks. 

1.1. Use in Education 

PIRATE needs to cater for a relatively large 
student body with little or no prior experience in 
observational astronomy, and with a very limited 
amount of remote supervision by a tutor. Grant- 
ing a sufficiently high level of access to the facility, 
such as the ability to initiate the opening and clos- 
ing of the dome, is considered essential for the stu- 
dent learning experience. The adopted software 
solution, a combination of the commercial prod- 
uct ACF0 (jPennv 2011) and the propri etary soft- 
ware described in iLucas fc Kofbl (l2010h . achieves 
the desired student access privileges without com- 
promising the safe operation of PIRATE. After 
one tutor-led induction evening conducted via au- 
dio communication and web interfaces, the stu- 
dent groups could operate PIRATE successfully 
on their own. A (remote) night-duty astronomer 
was on call but only occasionally contacted. 

PIRATE was deployed for OU undergraduate 
students for the first time in spring 2010 in a 10 
week project which is part of the third level (third 
year) OU module S382 Astrophysics. A total of 
about 30 students formed three groups with alter- 
nating access to PIRATE, for a total of 40 observ- 
ing nights. Each group sele cted a suitable target 
source from the catalogue bv lNorton et al. I ( 2007 ) 
of periodic variables found by Super WASP and co- 
incident with a ROSAT source. The groups then 
built up a long-term light curve of their target, 
such as the example shown in Fig. [TJ Individual 
observing sessions were staffed by small observer 
teams of 2-4 students who kept in audio and text 
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contact throughout the night. The collaboration 
of different observer teams in the larger groups on 
the same target source ensured the emergence of 
a useable data base to which all group members 
could develop a sense of ownership, even when 
occasional nights were clouded out and the ob- 
server teams on duty may not have succeeded in 
obtaining data themselves. A fuller account of the 
challenges and solutions for the PIRATE te aching 
project can be found in Kolb et al. ( 2010f) . 
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Fig. 1. — The light curve created from 206 I band im- 
ages of 1SWASPJ163420.90+424433.4. The data were 
taken to fulfill project objectives for the OU's under- 
graduate course, S382, by OU undergraduate students. 



1.2. Use in Research 

In its research role, PIRATE performs follow-up 
photometry of SuperWASP Transiting Extrasolar 
Planet (TEP) candidates from t he SuperWASP 
(Wid e Angle Search for Planets) (jPollacco et al 



20061 ) survey, and discovers extragalactic novae. 

Wide-field transit surveys produce light curves 
of millions of objects. Transit search algorithms 

are app lied to 
([20051 )) light 



(jCollier Cameron et al. | |2 006) 



trend-filtered (e.g. iTamuz et al 



curves to select TEP candidates from the millions 
of surveyed objects. Most of these candidates 
will not be TEPs, but will be 'astrophysical false 
positives' (AFPs), which are estimated to out- 



2 



number TEPs in a transit s earch survey by at 



least an order of magnitude (|0' Donovan et al 
2006). High priority candidates are followed up 



with radial velocity measurements to confirm the 
mass and therefore planetary nature of the transit- 
ing object, but this requires scarce large-telescope 
spectroscopy. 

To winnow TEP candidates from wide-field sur- 
veys like SuperWASP, small-to-medium sized tele- 
scopes provide higher precision , higher spacial res- 
olution follow-up photometry ( Haswell 2010f) . In 
this role, PIRATE can act as a link in the planet- 
finding chain, reducing the amount of large tele- 
scope time spent on false-positives. 

PIRATE has al so detected novae in M31 (e.g. 



M31N 2010-Ola) (|Burwitz et al. I l2010h . as part 



of an XMM-Newton large program for studying 
the source population in the Andromeda Galaxy, 
M31, in which a large number of supersoft X- 
ray sources were identified as counterparts of op- 
tical novae ( Pietsch 20081) . This requires the long 



term monitoring of M31 in X-rays with XMM- 
Newton, Chandra, and occasionally SWIFT. PI- 
RATE forms part of a wider network of optical 
telescopes that routinely monitor M3I to discover 
optical novae and track their light curves. Op- 
tical spectra of the freshly discovered novae are 
also obtained with larger telescopes, such as the 
Hobby-Eberly Telescope, USA, and the BTA, Rus- 
sia, to classify the novae. With the help of the in- 
formation extracted from these X-ray and optical 
data we are beginning to obtain a better picture of 
the physical processes that take place in nova out- 
bursts. This optical observing campaign grew out 
of the XMM- Newt on/Chandra M 3I nova monitor- 
ing collaboration (jPietsch 1 12010). 

In the following sections we describe the hard- 
ware and software used in the PIRATE facility, as 
well as detailing commissioning issues, data reduc- 
tion techniques, and future improvements. 

2. Hardware and Software 

PIRATE Mark 1 is a remote telescope facility 
consisting of a 14" Schmidt-Cassegrain reflector 
mounted on a Paramount MeH German Equato- 
rial Mount (GEM) (Fig. [5]). Its main imager is an 
SBIG STL-1001E, which houses the Kodak KAF- 



1001E CCD, a 1024 x 1024 pixel front-illuminated 
detector with a pixel size of 24/Lt, and a quan- 
tum efficiency of 0.4, 0.55, and 0.65 for wave- 
lengths 450, 550, 650nm respectiveljo In combi- 
nation with the 3.91m focal length of the Celestron 
C14, this produces a plate scale of 1.21" pixel -1 , 
and a field of view of 22' x 22'. An Optec TCF 
micro-focuser is located between the photometer 
enclosure and optical assembly, along with an 8 
position filter wheel, containing Johnson-Cousins 
standard BVRI broadband filters and 3 narrow- 
band filters (Ha, SB, Oil). PIRATE has a Cele- 
stron 80mm refractor with SBIG ST402 ME cam- 
era for auto-guiding. The telescope is housed 
within a Baader Planetariurr0 3.5m diameter All- 
Sky clamshell dome (Fig. [3]) with built in weather 
systems for automatic shut-down in adverse condi- 
tions. On-site weather is monitored by a weather- 
station and Boltwood cloud sensor (to the side of 
the dome), and both internal and external dome- 
mounted weather systems (rain, internal humid- 
ity, internal temperature). The dome has its 
own firmware which interfaces with the propri- 
etary dome driver, as well as with the Boltwood 
cloud sensor, allowing for shutdown conditions to 
be communicated directly to the dome, bypassing 
the control PCs. Four D-Link webcams (3 inter- 
nal, 1 external) provide live video and audio feeds, 
as well as IR beams for night- viewing. 




Fig. 2.— The C14 as seen through the D-LINK DCS- 
5230 IR webcam 
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the dynamic range of <1%, and <2% at the bot- 
tom end. We measured the gain using the fol- 
lowing relationship: a ADU = \(Nadu) (IHowell 




Fig. 3. — The 3.5m Baader Planetarium All Sky 
dome. The control room is directly below. 



2.1. Camera Characteristics 

To characterise the CCD we determined the 
gain, linearity, read noise and the dark current 
(which was measured at a range of temperatures) . 
The gain and the linearity measurements were 
made using the same set of dome fiats for each. 
A set of dome flats of increasing exposure time 
were taken in order to measure the median counts 
of a 100 x 100 pixel subframe in the middle of each 
image. For the short exposure times (of order a 
few seconds), we attempted to generate a shutter 
correc tion map using the methodology of IZissell ' 



(|2000l ). anticipating that we might see position- 
dependent corrections to the exposure time (due 
to the shutter travel) of order 10" 3 s. However, 
this produced a null result, and further investi- 
gation into the shutter mechanism of the STL- 
1001E confirms a rotating 'shutter wheel' that 
should be devoid of shutter travel effects. We 
therefore do not apply any shutter correction to 
the short exposure flat fields. We measured the 
bias level from contemporaneous bias frames and 
subsequently subtracted the mean pedestal level 
of 107 ADU from each flat. The subframe was 
chosen to be the centre of the vignetting function, 
where the image is at its flattest. To assess the 
linearity of the CCD, we plotted median subframe 
counts against exposure time (Fig. |4] (a)), and 
fit a linear trend to the same ADU range. The 
residuals of this fit can be seen in Fig. 2] (b). We 
note a deviation from linearity at the top end of 



2000), where a adu is the standard deviation of 



the sub-frame counts, and (Nadu) is the mean 
recorded counts in the subframe. To derive this 
expression, one assumes the statistical relationship 
°e- = VN e - holds, where a e ~ and N e - are the 
uncertainty in and the number of recorded photo- 
electrons respectively, which only holds when the 
process is governed by Poisson statistics. To de- 
termine g we therefore plot <J 2 ADU against (Nadu) 
(see Fig. ffl(c)). We needed sufficient photoelec- 
trons for photon-counting (Poisson) noise to be 
dominant, and so we fitted to the range 20k ADU 
< (Nadu) < 40k ADU. We chose an upper limit 
of 40k ADU to be well clear of the digital count- 
ing limit of 65535 ADU (16-bit). We measured a 
value for the gain of 1.62 ± 0.01 e-ADU" 1 . The 
manufacturers quote a full- well capacity of 150k 
e _ , which, given the determined gain, would oc- 
cur at a theoretical ADU value of ~ 92.5k ADU, 
an unobtainable value due to the digital counting 
limit. We note that falling short of the electron 
full well capacity is probably responsible for our 
excellent linearity of response at the upper end of 
our digital counting range. 

To measure the read-out noise, we took 30 
bias frames at a temperature of —20°, and me- 
dian combined them. The median-combined bias 
frame (which we assume to have negligible read 
noise) was then subtracted from each of the 30 
bias frames, and the standard deviation of each of 
these difference images measured, we aver aged this 



figure to achieve the read noise in ADU ( Howell 



2000) . We measure then (using the previously de- 
termined gain value) a read-out noise of 10.9 ± 1.3 
e _ . To check for consistency, we also split the 
bias frames into pairs, subtracting one from the 
other within each pair, so that d = b\ — &2> where 
b\ & 62 are the bias frames in each pair, and 
^read = &d/v2, yielding a read- noise measure- 
ment for each pair. These values were consistent 
with the previous read-noise measurement. 

Dark current was assessed for a range of tem- 
peratures (from 15°C through to -25°C in in- 
crements of 5°C) by median-combining 10 darks 
at each temperature, 10 biases at each tempera- 
ture, subtracting the bias and dividing the resid- 
ual counts by the exposure time (30s for each 
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dark frame). The measured values are listed be- 
low in Table [T] We measure a dark current of 
0.02 ± 0.02 e-s" 1 and 0.08 ± 0.02 e's" 1 at -20°C 
and -10°C, the typical operating temperatures for 
winter and summer respectively. We also find a 
value of 0.51 ± 0.03 e^- 1 at 0°C, c.f. the 1 c^ 1 
at 0°C manufacturer's specification. 

2.2. The Mount 

The Paramount ME GEM is able to deliver an 
all-sky pointing accuracy of around 2'. We em- 
ploy the software 'TPointH to model system flex- 
ure to achieve this value, in which we use 30 model 
points that capture the discrepancy between in- 
tended and actual pointing for a variety of tele- 
scope positions. We had previously employed a 
~150 point model, which was abandoned due to 
the addition of new hardware. We noted that the 
subsequent 30 point model replacement provides 
sufficient pointing accuracy to enable the control 
program to plate-solve and correct for pointing in- 
accuracy with only a small contribution to inter- 
frame overhead, and so have settled at 30 points 
for convenience. We also have the ability to iden- 
tify and remove periodic error in the sidereal track- 
ing from within the mount control software, by 
measuring the correcting impulses provided by the 
autoguider. The periodic error is caused by ir- 
regularities in the rotation rate of the worm gear 
through one cycle, and so presents itself entirely 
in RA. We fitted a 5 th order polynomial with a 
peak-to-peak amplitude of 4.8", and subsequently 
saw improvements in the circularity of each star's 
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T (°C) 


D. C. (e" s- 1 ) 


a d.c. (e- ADIT 1 ) 


-25 


0.01 


0.02 


-20 


0.02 


0.02 


-15 


0.05 


0.02 


-10 


0.08 


0.02 


-5 


0.27 


0.03 





0.51 


0.03 


5 


1.13 


0.03 


10 


2.25 


0.03 


15 


5.20 


0.04 



Table 1: Measured dark current for a range of chip 
temperatures 




0123456 012 3 456 

Exposure time (s) Median Subframe Counts (10 4 ADU) 




12 3 4 5 6 -30 -20 -10 _ 10 

Exposure time (s) CCD Temperature ( C) 



Fig. 4. — (a) Top left: Counts as a function of 
exposure time, with a linear fit through the expected 
linear region 20k - 40k ADU. (b) Bottom left: the 
residuals of the above fit. (c) Top right: the photon 
transfer curve, with linear fit to the same 20k - 40k 
ADU region. The deviation from linearity at ~0 - 
10k ADU is due to the transition from the shot noise 
to the read noise regime, (d) Bottom right: dark 
current as a function of chip temperature 
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point spread function. As this is a German Equa- 
torial Mount, the OTA must 'flip' from one side 
of the pier to the other when the tracked target 
crosses the meridian. This results in the observed 
field being rotated 180° with respect to the optics 
and camera. 

2.3. Software 

The hardware is controlled by a collection 
of Windows-based software (Maxim DL@ for 
ima ging, The Skyd*"! fo r mou nt control, Focus- 
MaaPlf Weber fc Brady I l200lh for control of the 



micro-focuser) which are each controlled by ACPi 12 !. 
an observatory control suite that automates the 
command of each of the component software pack- 
ages through the ASCOM-standard driver layer. 
These programs are run on a 2.13GHz CPU Win- 
dows XP PC with 2Gb RAM. The site receives 
its internet connection via satellite; connection 
speeds are good, with typical upload speeds from 
the OAM to the UK ~ 250kBs~ 1 , and operating 
latencies of around 60-100ms. We use Dimension 
4 from Thinking Man Software 

E 

to synchronise 

the control PC's clock to UTC via an NTP server 
every 10 seconds. 

For a given observation or group of observations 
in a night, ACP will first typically initiate an aut- 
ofocus via FocusMax, which searches for mag 4-7 
stars within a 2° x 2° region to focus with. Once fo- 
cused, ACP starts a 4s 'pointing exposure', which 
it subsequ ently plate-solves vi a PinPoinO using 



the GSC (jLasker et al. I Il996) to determine the 



pointing error. A corrective slew is applied, and 
the image sequence commences. Each full expo- 
sure, once read-out, is plate-solved once more to 
check the pointing. We apply a maximum pointing 
error of 3", forcing a pointing update if a given im- 
age is determined to be misaligned by more than 
3". The benefits of doing so can be seen in Fig- 
ure [5J where the centroid deviation is significantly 
reduced. This reduces the contribution to photo- 
metric noise from flat fielding errors; more on this 
in section [331 



9 http: / /www. cyanogen, com 
10 http:/ /www. bisque. com 

11 http:/ /users. bsdwebsolutions.com/ larryweber/ 
12 http:/ /acp. dc3.com 

13 http: / /www. thinkman.com/dimension4 
14 http:/ /pinpoint. dc3.com/ 



The software tasks performed between expo- 
sures in a time-series run are: 

1. The previous image is read out and plate- 
solved, from which the astrometric residual, 
pointing error (magnitude and direction), 
mean FWHM, true focal length, true image 
centre, and camera sky position angle are 
calculated. 

2. The centering is tested: if within the max- 
imum pointing tolerance, no slew update is 
performed. If outside the maximum error 
tolerance, the auto-guiding is temporarily 
stopped, and a corrective slew is performed. 

3. A guide star is selected from the guider im- 
age, with the star providing a signal-to-noise 
> 3 for the shortest exposure interval chosen. 
The exposure is started. 

We typically see astrometric residuals from the 
plate-solve of ~0.2". The total inter-frame over- 
head between exposure end and exposure start is 
~ 15s. Each frame takes ~ 2.5s to read-out and 
download, meaning the overhead is dominated by 
the time taken to find the plate solution. The 
15s quoted refers to the case where no subsequent 
pointing update is required after the plate-solve. 
When a pointing update is required (pointing er- 
ror >3"), l-2s is added to the overhead time. 

ACP also provides a web control interface, 
which made it an obvious choice for implementa- 
tion in a distance-learning, group-based education 
scenario. When used by researchers, direct access 
to the PIRATE control PC is achieved by remote 
access of the control PC's desktop. Students, 
however, are not given direct access to the main 
control PC, and instead interact with the hard- 
ware through a proprietary server/client set-up 
that limits their control. The server-side appli- 
cation ('Switch Server') provides access to relay 
switches which are used to turn components on 
individually, and also has a suite of extra ca- 
pabilities. These include: preventing the dome 
from being opened in daylight hours; closing the 
dome in case of excessively high humidity or strong 
winds; automatically cooling the two cameras in 
a monitored, staggered sequence; and checking 
the availability of the local internet connection, 
forcing a dome shut-down in case of lost internet 



6 




-10 -5 5 10 15 20 
Centroid deviation from mean in x (pixels) 

Fig. 5. — Position variations during a night with no 
specified maximum pointing error (black circles) and a 
night with maximum pointing error set to 3" (crosses). 



connection. Further information on the opera- 
tion of PIRAT E's Sw itch Server can be found in 
(1201 0l). 
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3. Data Acquisition, Reduction & Com- 
missioning Issues 

3.1. Observing Strategy 

On a given night, the user cools the main imag- 
ing camera down automatically using the Switch 
Server to a target temperature of ~ 40° C below 
the current ambient temperature. Consequently 
the operating temperature of the CCD varies sea- 
sonally, usually between -10°C in the Summer and 
-20°C in the Winter. This ensures the camera is 
always operating at the lowest available temper- 
ature. Once the target temperature is achieved, 
ACP is started and the dome opened around sun- 
set. ACP's automatic sky flat procedure is then 
initiated. This procedure automatically scales ex- 
posure times to achieve a target peak count of ~ 
33k ADU in each filter's flat fields. During the flat 
procedure, sidereal tracking of the mount does not 
occur. The flat fields are dithered so that any stel- 
lar sources present in the sky flats can be removed 



upon combination of each flat into the master flat 
field. Once the flat fields are obtained, the dome 
is then shut whilst waiting for the end of Nautical 
twilight (ACP allows observations through Astro- 
nomical twilight). During this time, the user takes 
calibration frames. 50 bias frames are taken first, 
to assist in completely flushing the chip of any 
residual charge from the flat fields. Then dark 
frames are taken in batches of 10 for each antic- 
ipated exposure time to be used for the science 
frames throughout the night. Typically these ex- 
posure times are 30, 45, 60, 90, and 120s. The 
dome is then reopened in time for the end of nau- 
tical twilight, and the observing schedule for the 
night commences. 

3.2. Photometric Reduction, Performance, 
and Light Curve Generation 

The data is available in real-time via direct 
FTP to the control PC, and users with direct ac- 
cess to the control PC may also inspect each image 
'live' after being read-out from the camera. Each 
morning the previous night's science and calibra- 
tion frames are then transferred from the OAM to 
a local archive at the OU; the user may obtain the 
frames the following day. Reduction (in the re- 
search case) is performed by a custom IDL script 
that operates interactively, performing file renam- 
ing and sorting tasks and determining the pres- 
ence of calibration frames, subsequently compiling 
IRAfB scripts that ma ke use of the DAOPHOT 



package (Ste tson 1 11987t h The DAOPHOT out- 



put is then fed back into an IDL script that per- 
forms differential photometry and generates the fi- 
nal light curves for each object in the frame. When 
used for teaching, the students do not make use of 
such scripts, and must process their frames 'man- 
ually' using Maxim DL. 

The initial IDL script first compiles a list of 
available light frames from the previous night, and 
then iteratively displays each one in turn, asking 
the user to reject or accept each frame. This allows 
for the rejection of frames with defects, such as a 
lack of useable point sources in the frame, trails, 
or excessive pointing error. The user then accepts 



15 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 
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or rejects a final list of rejected frames, producing 
a master file list on acceptance. Lists are gener- 
ated for each calibration frame type; discriminat- 
ing flat fields by filter. If flat fields matching the 
chosen filter of the light frames do not exist, no 
flat fielding is performed. An IRAF script is gen- 
erated for the task CCDPROC, which performs 
standard CCD calibrations and then executes the 
requisite post-calibration rotation commands for 
pre-flip light frames. 

The IDL routine first proposes a master frame 
according to a simple statistical test performed 
on all frames. The master frame proposed will 
be the frame that exhibits the highest standard 
deviation from the quartile of frames that have 
the lowest mean values. The mean test ensures 
a low sky background (high sky backgrounds are 
often indicative of Cirrus scattering) , whilst select- 
ing a frame with high standard deviation ensures 
many strong point source fluxes well above the 
background. The user can over-ride the proposed 
frame selection and select their own frame as the 
master frame. An alternative method would be 
to determine the frame with the minimum aver- 
age FWHM, and designate this the master frame 
due to the superior seeing. This is not currently 
implemented. 

Once defined, we use DAOFIND on the master 
frame (making use of a user-given average frame 
FWHM) to find all sources above a 6a threshold. 
An initial PHOT task run on the master frame 
determines the master coordinate list, with only 
sources that do not suffer a PHOT task failure 
flag going into the master coordinate list. A lin- 
ear transformation is then deduced between the 
master frame coordinates and those of all other 
frames, and this is used to generate the final co- 
ordinate lists for each frame. Final photometry is 
standardly performed with an aperture 3 x the in- 
put (master frame) FWHM, though the aperture 
size can be changed if the data warrants (e.g. in 
the case of close companions to the target that 
lie within or partially contaminate the measuring 
aperture) . 

Once the instrumental photometry is obtained, 
it is fed into another IDL script which places the 
whole night's photometry in memory as a data 
cube. Mid-exposure time in UTC then is con- 
verted to HJD for all frames. We note the prefer- 
ence of using BJD over HJD, which can differ by 



as much as 4s ( Eastman et al. 2010f) . as well as 
calculating BJD from TT instead of UTC, which 
currently differ by 66.184s as of July 2011. We will 
be incorporating these changes into IDL scripts in 
the near future. The night's data is then investi- 
gated for PHOT task failures. Light curves that 
suffer a photometry failure in one or more frames 
are rejected from the cube. In case of bad frames 
still present in the data (which would lead to an 
excessive number of light curves being cut), we in- 
troduce a parameter (to be set by the user) defined 
as: Ngtar/a, where a frame is removed if it con- 
tains a fraction of photometry failures > N star /a. 
Here N star is the number of stars in the frame, and 
a is a free parameter. We find a value of a = 10 to 
usually be sufficient. The cube is then sorted by 
magnitude, and an iterative ensemble compilation 
procedure begins . We employ a simi lar method- 
ology to that of Burke et al. ( 20061 ). where for 
each star, all the corresponding comparison light 
curves are sorted according to a light curve 'fig- 
ure of merit'. Each comparison light curve is nor- 
malised to unity, and these light curves' RMS (i.e. 
the standard deviation about the mean) is used as 
the figure of merit. Starting with the best rated 
light curve (i.e. lowest RMS) for a given star, 
each subsequent light curve in the list is iteratively 
included in the ensemble, using inverse variance 
weights, where the weights derive from the IRAF 
photometric uncertainties. Note that the ensem- 
ble never contains the flux from the target star, 
as the target star is excluded from the list of light 
curves. The following refers to the production of 
the final light curve for one star only. The en- 
semble flux, Ej, for a given frame j, containing 
comparison stars i = 1...N is compiled as thus: 



E, 



Where Fji are the stars' individual fluxes, and the 
individual weights are given by: 



1 



The uncertainty in the ensemble is given by the 
standard deviation of the weighted average: 



' N 

E" 



-1/2 
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With each new star added to the ensemble, the 
differential target star light curve d* = (d\ . . . djv), 
where d^j = = 1 . . .N, is calculated from 

the recorded flux of the target star, F* j, and cal- 
culated ensemble flux, Ej. The light curve d* is 
then normalised by dividing through by its me- 
dian value, and the RMS is computed. The fi- 
nal light curve, d„ji n makes use of the ensemble 
that produces the lowest RMS. Each star there- 
fore has its own unique comparison ensemble. In 
Figure O we display the results of this method 
for the stars in the WASP- 12 field on two con- 
secutive nights with approximately 30% & 39% 
moon illumination, at a separation of ~ 142° 
and ~ 132° fro m the moon respectiv ely, compar- 
ing USNO-B1.0 (|Monet. et al. ll2003l) magnitudes 
with the RMS precision of PIRATE R band light 
curves from 60s exposures. The observation runs 
of 22/11/2009 and 23/11/2009 contain light curves 
of 759 and 730 objects, each consisting of 359 and 
321 exposures respectively. We see from this plot 
that sub-percent precision is available for magni- 
tudes ~13.5 and less with an operating cadence 
of 76s (including the aforementioned overheads). 
We include for comparison the contribution to the 
uncertainty from photon noise alone. 
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Fig. 6. — Light curve RMS deviation as a function of 
USNO-B1.0 R magnitude for the nights of 22/11/2009 
and 23/11/2009, 60s exposures, WASP-12 field 
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Fig. 7. — WASP-12b transit observation (R band) 
from 23/11/2009. Residuals from the model fit are 
shown offset by 0.05. 



As an example, this particular method of light 
curve compilation was used in the production 
of light curves of the transit of WASP-12b, for 
which 6 nights of data were captured during the 
months of October-November 2009, the details of 
which can be found in Haswell et al. (in prepara- 
tion). We show the light curve from the night 
of 23/11/2009 in Figure which has a post- 
fit residual RMS of ^2.6 mmag. This data was 
used to confirm the published ephemeris for an 
HST visit investigating t he exosphere of WASP- 



12b dFossati et al. 1 120101 ) . The displayed model 



light curve is the published model light curve (fit 
to the data from the disc o very p aper) using the pa- 



rpa. 

rameters of iHebb et al. I (120091) . and analy tic light 



curve solutions of iMandel fc Agol I d2002l). using 
the limb-darkening coefficients of IClaret I (|2000l) . 

We find that the entire light curve's RMS is 
an unsuitable figure of merit for determining how 
accurately a light curve has captured an object's 
inherent variability in certain cases. Substantial 
out of transit coverage of WASP-12 (out to phase 
1.4 (not shown)) meant that the RMS is an ac- 
ceptable figure of merit for the observations shown 
in Fig. [7J Obviously, if one has prior knowledge 
of the transit ephemeris, as with these WASP-12 
observations, solely the out-of-transit light curve 
should be used to calculate the figure of merit. 
In the case of the follow-up of candidate plan- 
ets, or for other purposes (e.g. cluster surveys), 
the ephemeris might not be well constrained. The 
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method's success is therefore subject to the ampli- 
tude, timing, and character of the star's intrinsic 
variability. 

When tackling variable stars with ampli- 
tudes over and above ~10%, the procedure 
will often select very faint, non-varying stars 
for inclusion within the ensemble, so long as 
their light curve RMS is less than the intrin- 
sic variability of the target star. We showed 
in Fig. Q] the light curve of a larger amplitude 
(~20%) vari able, 1SWASPJ163420.90+424433.4 



([Norton et al. I 120071 ). taken on 07/06/2010 by 



OU undergraduate students as part of their prac- 
tical module. The data were reduced separately 
for inclusion here, with 4 comparison stars cho- 
sen by hand for combination through inverse- 
variance weighting into a comparison ensemble, 
as the RMS-minimisation method does not work 
for this target. We also find the procedure inade- 
quate in some lower amplitude situations, such as 
a typical hot jupiter transit with insufficient out- 
of-transit coverage. For example, a light curve 
that captures 1 hour in-transit, egress, and 1 hour 
out-of-transit approximates a step function; and 
is bimodal in d* , a scenario for which the RMS is 
a poor figure of merit. 

To combat this, we are developing an exten- 
sion of the technique that breaks each comparison 
light curve into multiple segments, and performs a 
linear regression in each segment, producing a se- 
lection of relevant statistics: segment RMS, resid- 
ual segment RMS (after subtracting the fit), seg- 
ment median value. By careful analysis of each 
light curve's segment statistics, as well the statis- 
tics of the segments as a whole (e.g. range or 
standard deviation of the segment medians), it 
will be possible to identify variable stars for ex- 
clusion from a comparison ensemble, whilst simul- 
taneously allowing for any variable star to have an 
optimal comparison ensemble automatically deter- 
mined. This technique will facilitate serendipitous 
PIRATE discoveries of previously unknown vari- 
able stars. 

3.3. Flat Field Inaccuracy and Optical Re- 
sponse Calibration 

The systematic errors of fiat fields as optical 
response calibrators are often overlooked in the 
application of differential photometry, as modern 
auto-guiding systems are sufficiently good at lock- 



ing the stellar point-spread functions to fixed po- 
sitions on the chip for the duration of a time se- 
ries. As the observer is only interested in time- 
variability, systematic offsets in the magnitude 
zero-points of any two stars in the field can be 
ignored. PIRATE's GEM executes a 'pier-flip' 
as the tracked target crosses the Meridian. This 
moves the OTA to the other side of the pier; and 
inverts the image of the stellar field with respect to 
pre-fiip frames. This operation effectively moves 
each stellar point spread function to an entirely 
different part of the focal plane. To maintain con- 
tinuity in the flux ratio between two stars across 
the pier flip, the vignetting must be determined 
perfectly. In reality, most flat fields suffer from 
inaccuracies at the 10~ 2 level (|Manfroid I [l995l) . 



We therefore expect and indeed see varying flux 
ratios between two objects in the field across the 
pier-flip. We term this effect a 'Light Curve Dis- 
continuity' (LCD). Several mechanisms contribute 
to the flat field systematic errors: 



Uneven illumination from the source. This 
can occur in sky fl ats, even when they image 
the sk y null point ([Chromev fc Hasselbacher 
19961 ) as PIRATE's automatic procedure 
does. It will also occur in dome flats. 



Physical changes in the hardware between 
taking flat fields and science images that 
modify the response function of the system. 
In the case of PIRATE Mark 1, it is sus- 
pected that primary mirror 'flop', which oc- 
curs despite the presence of mirror locks, 
and/or flexure in the micro- focuser assem- 
bly changes the vignetting function. It may 
also lead to the misalignment of the focal 
plane with the CCD surface, meaning con- 
sistent focus across the field of view is not 
achievable. In this scenario, the FWHM of 
the stellar profiles will be a function of chip 
position, and the measuring apertures used 
for photometry would see different levels of 
flux leakage due to the changed stellar profile 
after a pier-flip. 

Scattered light in the optical tube assembly. 
Off-axis, unfocused stray light due to insuffi- 
cient baffling can contribute to the flat field 
exposure. Scattered light is additive. If in- 
cluded in the flat field, this additive light 
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is incorrectly used in the multiplicative flat 
field calibration. 

We have been unable to develop a wholly effec- 
tive procedure for self-consistently calibrating the 
LCD effect without investing much of the night 
in taking calibration observations. We thoroughly 
explored these possible approaches: 

We used data from the night of 23/07/2009, 
during which the moon had 3.2% illumination, so 
we might expect sky background gradients to be 
limited. The data were processed in the usual 
manner, and investigated separately as two groups 
of 'pre-fiip' and 'post-flip' frames. The frames in 
each group were median-combined, and a 6 th order 
Legendre polynomial was fit to the sky background 
in each of the two (pre and post) resultant frames 
using the IRAF task 'IMSURFIT'. Note that the 
pre-flip frames were rotated 180°. We denote 
the background fits by A pre (x,y) and A post (x, y), 
where x and y are image co-ordinates correspond- 
ing to positions on the sky (not pixel co-ordinates, 
due to the aforementioned rotation of the pre-flip 
frames). A simple ratio of the pre and post-flip 
background fits reveals any discrepancy in the sky 
background for a given star position. We denote 



this ratio as A map {x,y) 



A pr<1 (x,y) 
t{x,v) ' 



We show 

Ama P in Fig. [U We note that it has structure pre- 
dominantly in the x direction, and displays a peak- 
to-peak variation of ~ 7%. The (non-differential) 
light curves of each of the N stars (here N = 517), 
which we denote by Fi(t), where i = 1...N were 
median-combined to create an approximation (to 
first order) of the sky transparency, F. We assume 
each star remains at a fixed position in the image 
(xi,yi) for the duration of the observing run. For 
all of the light curves (including the transparency 
function), the fluxes were averaged over time for 
pre-flip frames and also for post-flip frames, pro- 
ducing a single pre-flip and post-flip flux value 
for each light curve. We therefore use F pre ,i and 
F pos t,i to refer to these time-averaged values. We 
define the LCD of a star to be: 



AR = 



»/# P 



We now have the ability to assess the applica- 
bility of the background 'map'; if the shape of the 
sky background (after flat-fielding) provides an es- 
timation of any residual error in the flat-fielding 



process, then it should correlate well with the ob- 
served LCDs. To compare the map with the LCDs, 
we fit a 2 nd order polynomial through least squares 
regression to the irregularly gridded AFj values. 
The result of this fit can be seen in Fig. |H1 As can 
be seen, far from correlating well, the pattern of 
the LCDs also exhibits greatest deviation in the 
x direction, but with opposing orientation. The 
structure of the map in Fig. |8] suggests a fixed 
and constant light source co-moving with the op- 
tical tube assembly, inducing the same background 
structure in both pre and post-flip frames, which 
is then amplified when the pre-flip frames are ro- 
tated. One would expect any transformation of 
the vignetting function to be recorded in the sky 
background, as this gives a continuous indication 
of the vignetting function. However, the presence 
of any scattered light in either the science frames 
or flat fields renders the true vignetting function 
for each side the pier flip unrecoverable. 

The scattered light structure present in the sci- 
ence frames could also have been present in the 
flat field, and thus the LCDs introduced to the 
data through the flat fielding process alone, but 
the un-flat fielded data shows the same LCD struc- 
ture seen in Fig. [S] Given the dark night (one 
night after a new moon) for this data set, it is 
apparent that scattered (additive) light is likely 
always present, and greatly reduces the accuracy 
of standard flat-fielding and sky-fiat procedures; 
the GEM simply highlights this inaccuracy. 

The data of 23/07/2009 exhibit strong struc- 
ture in the LCD map, and is hence a good demon- 
stration of the problem. A more typical structure 
can be seen in Figs. [TUlandrm from the night of 
23/11/2009, WASP-12 field. The polynomial fit 
to the LCDs exhibits a smaller peak-to-peak am- 
plitude of ~2%, and there is significant scatter in 
the residuals of this fit to AFi . 

We conclude that background fitting, sky flats, 
and even flat fielding are all ineffective at negating 
the percent-level LCD effect seen in PIRATE Mk. 
1. Instead, there are two routes to achieving cor- 
rect photometric calibration across the pier-flip. 

The first of these involves creating a pho tomet 
ric sup e r-flat. This is d i scusse d in detail in Boyle 
" (l2008lL iManfroidl (Il995h . 



(2007), Grauer et al. 



Regnault et al. I (|2009), and ISehnan I ([20041 ). At 



its most simple, this involves observing a set of 
standard stars at different positions in the field of 
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Fig. 8. — Sky background map produced by fitting to 
the sky background for all pre-flip frames and all post- 
flip frames, then taking the ratio of these two fits. The 
structure of this map suggests there may have been a 
fixed scattered light source with respect to the OTA. 
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view, and determining the position-dependent re- 
sponse to the standards in order to build the large- 
scale, low-frequency optical response of the system 
into the final flat field. Sky flats or dome flats 
are observed in tandem, with polynomial fits to 
their large scale variation made in order to 'flatten' 
them. This retains the small scale, high frequency 
intra-pixel response, which can then be combined 
with the standard star-determined low frequency 
response function to create a 'photometric super- 
flat'. In constructing the large-scale component, 
it is preferable to make multiple dithered obser- 
vations of a cluster containing many standards 
that span the field of view in order to reduce the 
number of frames needed. As we expect the re- 
sponse to be different from one side of the pier to 
another, this involves observing the same cluster 
twice, once in the eastern sky, and again later in 
the western sky. 

The second method involves lending a degree 
of freedom to the normalisation of each light 
curve section (pre & post-flip), shifting the flux 
levels up and down to minimize the \ 2 °f a 
model fit. This method has its limitations in the 
extent of fore-knowledge required to work suc- 
cessfully. In the case of SuperWASP follow-up, 
the Markov-Chain Mon t e Car lo fitting routine of 
Collier Cameron et al. ] f)2007h can take the light 
curves from each side of the pier as separate input 
light curves and include the normalisation factor 
as a free parameter in the fitting procedure. How- 
ever, this algorithm works with 'prior knowledge' 
of the light curve, in that its initial parameters 
from which the MCMC routine iterates are taken 
from the SuperWASP light curves. This prior ex- 
pectation of the light curve's model parameters 
(including, crucially, the transit ephemeris) as- 
sists with the cross-flip normalisation. If there is 
no prior knowledge of the light curve, adjusting 
flux levels is almost certainly perilous. Take, for 
example, the worst case scenario of a pier-flip oc- 
curring mid- way through ingress or egress. In such 
a scenario it would prove difficult to determine if 
the flux deficit across the pier-flip is instrumental 
or astrophysical. Great caution is called for in 
interpreting such light curves. 



4. Summary and discussion 

PIRATE is a remotely-operable telescope facil- 
ity built primarily from readily available 'off-the- 
shelf components, used in both education and re- 
search. In its education role, it has successfully 
supported small groups of simultaneous student 
users in their efforts to observe and classify vari- 
able stars in the SuperWASP archive. We showed 
that relatively inexpensive equipment can play a 
useful role in modern astrophysical research, in- 
cluding high-precision time-series photometry re- 
quired for transiting Jupiter-size exoplanets. The 
PIRATE system employs a GEM which provides 
exceptional stability and all-sky pointing accuracy 
for its cost. The pier-flip introduces LCDs which 
are difficult to correct for by calibration proce- 
dures alone; the best strategy is to treat pre- and 
post-flip light curves as separate observing runs. 

The OTA and camera combination described 
here awaits re-assembly in a neighbouring dome 
at the OAM. We note that the combination of a 
Paramount ME and the affordable Celestron C14 
can be found in many installations world-wide. 

The use of a remotely-operable telescope in the 
OU module S382 proved the feasibility of de- 
ploying such complex hardware in real-time to cre- 
ate an inspirational teaching tool for distance ed- 
ucation. The great enthusiasm of the students in- 
volved, and the quality of the acquired data and of 
the scientific reports generated at the end of the 
project demonstrate the success of the PIRATE 
teaching project. Use by level 2 (year 2) students 
is being implemented, and use by level 1 (year 1) 
OU students is currently under consideration. 

PIRATE has very recently undergone an OTA 
change to a Plane Wave Systems CDK17 f/6.8 Cor- 
rected Dall-Kirkham astrograph, which provides a 
bigger field of view and larger aperture (0.425m). 
Future work will involve linking ACP to the dome 
control via the ASCOM driver layer, allowing for 
full scripting of a night's events; adding a level 
of autonomy to the observing process. Use of 
ACP ScheduleJ^l with the current weather sta- 
tion hardware and switch-server intelligence would 
effectively allow PIRATE to run in a fully au- 
tonomous, robotic mode, which will be explored 



16 http: / /www3. open, ac.uk/study/undcrgraduate/course/s382. htm 
17 http:/ /scheduler. dc3.com/ 
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as an operation mode for research purposes. The 
real-time remote operation will remain the only 
mode of operation for teaching applications, how- 
ever, in order to make the learning process as 
interactive as possible. 

We thank the referee for their constructive 
and insightful comments, which improved the pa- 
per. We are grateful for the continued support of 
Baader Planetarium towards optimising the dome 
for the PIRATE facility, and the on-site support 
provided by staff at the OAM. SH acknowledges 
the support of an STFC studentship. 
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